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Abstract: We investigate possible signatures of a pre-infiationary anisotropic phase in two- 
point and three-point correlation functions of the curvature perturbation for high-momentum 
modes which exit the horizon well after isotropization. In this momentum regime, the early 
time dynamics admits a WKB description and the late time dynamics can be described in 
terms of a non-Bunch Davies vacuum state which encodes the information of initial anisotropy 
in the background spacetime. We compute the bi-spectrum for curvature perturbation in a 
canonical single-field action with and without higher derivative operators. We show that the 
bi-spectrum at late times, in either case, is enhanced for a flattened triangle configuration 
as well as a squeezed triangle configuration and compute the corresponding /jvl parameters. 
The angular dependence and the particular momentum dependence of the /jvl parameter 
appear as distinctive features of background anisotropy at early times. 



Keywords: Anisotropy, Power Spectrum, Bi-spectrum. 



Contents 

1. Introduction 1 

2. Background Equations of Motion 2 

3. Cosmological Perturbations: Small Wavelength Limit 4 

4. Computation of the 3-point function 8 

5. Conclusion 15 

6. Acknowledgements 17 

7. Appendix: Scalar Field in Anisotropic Background 17 



1. Introduction 

A deeper understanding of the inflationary scenario requires going beyond the power spec- 
trum and probing higher correlation functions, collectively referred to as non-Gaussianitiy 
[1]. Observable non-Gaussianity, in turn, requires a departure from the standard single-field 
inflation with a canonical action [2, 3]. Substantial progress has been made in understand- 
ing the enhancement in non-Gaussianity for several variants of the standard scenario, which 
involve having multiple scalar fields, non-canonical action for the scalar field, introducing 
higher derivative terms in the action or having a non-standard vacuum state (see [4] for re- 
views) . A common feature of all these models is that they have a homogeneous and isotropic 
background for the perturbations to evolve. 

In the present work, we investigate the possible signatures of a strong anisotropy in the early 
space-time metric in the perturbation spectrum and its non-Gaussianity. The effect of pri- 
mordial anisotropy in the power spectrum has been studied before by [9, 8] and most recently 
by [10]. We redo the computation here for completeness. A primordial anisotropy in the 
metric will be washed away by a period of inflation, hence it seems interesting to try to find 
signatures that do not get washed away. 

We focus on a family of axially symmetric Bianchi I background geometries which admit 
a WKB solution for the perturbations at early times for modes in the high-momentum regime. 
On matching the WKB solution with the solution at late times, we can describe the late time 
dynamics of the curvature perturbation in terms of a non-standard ground state (essentially 
an excited state on the BD vacuum). The possible enhancement of the non-Gaussianity for 
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non-standard (non Bunch-Davis) vacuum was pointed out years ago by Holman and Tolley 
[5] . By studying the three-point correlation function in the present scenario and deriving the 
relevant contribution of background anisotropy to the f^L parameter, we infer that a possible 
enhancement in the bi-spectrum may also occur in the squeezed triangle limit, in addition to 
the flattened triangle limit discussed in [5]. 

Although we are primarily interested in non-canonical vacuum states which arise as a result 
of spatial anisotropy at early times , our computation easily generalizes to any excited state 
obtained by a Bogoliubov transformation on the Bunch-Davis vacuum. 

The paper is organized as follows. In the next section, we discuss the anisotropic background 
geometry in which the inflaton evolves and study the classical equations of motion. Section 
3 deals with the study of cosmological perturbations in the high-momentum regime and the 
WKB solution. In section 4, we compute the bi-spectrum first for a canonical action and then 
in presence of higher derivative terms. 

While we were preparing this manuscript, we came across the work [11], which has some 
overlap with our work. 



2. Background Equations of Motion 

We consider a theory of Einstein gravity with a minimally coupled single scalar field given by 
the following action, 

S = \ j d 4 xV=9~R + j d 4 x^ (-^d^d u( f> - V(cf>)\ , (M 2 p = 1) (2.1) 

where the background metric is chosen to be an axially symmetric version of the Bianchi I 
metric: 

fa* = _ dt 2 + e 2p (dx l)2 + e 2/3 (dx a )2 {%2) 

with a = 2, 3. 

In contrast to the FRW case where one has a single Hubble constant, we have two Hubble 
constants, which we choose to define as follows: 

*=£±M, h = t^ (2 . 3) 
3 ^3 

The classical dynamics of the system specified by the action (2.1) constitutes a strongly 
anisotropic expansion at early times (parametrized by h) followed by eventual isotropization 
at a time-scale t « ti so = For t » ti so , the universe enters a phase of de Sitter expansion. 

Note that h, which, roughly speaking, is a measure of the rate of anisotropic expansion 
vanishes in the isotropic limit (p = $) so that we are left with a single Hubble constant. 
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In terms of H and h, the independent Einstein's equation and the equation of motion for 
the scalar field reduce to the following set of equations: 



H + 3H 2 = V{4>) 
3H 2 - h 2 = + V(4>) 
4> + 3H<p + V'((f)) = 

The time-evolution of h can easily be derived from the above equations, 

h(h + 3Hh) = 

In the anisotropic phase, h ^ 0, which leads to the equation of motion, 

h + 3Hh = 



(2.4) 
(2.5) 
(2.6) 

(2.7) 

(2.8) 



For a general V (</>), one can only obtain approximate solutions to the above system of 
equations. However, for a constant V, one can exactly solve the coupled differential equations 
for H, h and <p as follows, 



H = 



coth(vWt) = Hi coth(v / 3W) 



h = ±W- 



sinh(V3yt) 



(2.9) 



= 



In the above solution, the constants have been chosen such that the metric approaches a 
Kasner solution in the limit t — > + . The ± sign in the solution of h indicates two different 
branches in the solution space (distinguished, among other things, by their behavior in the 
Kasner limit). It turns out that only for the positive branch, one can impose initial con- 
ditions on the cosmological perturbations at early times via the usual WKB approximation 
[12]. Hence, in this note, we will focus exclusively on this class of backgrounds. 

Now, for a given non-trivial V((p), the slow-roll condition ((f) 0) will imply that at early 
times V(4>) is nearly constant with time, provided Hcj) 2 — > at early times. This condition is 
obeyed by all common inflaton potentials and hence the above solution (2.9) can be trusted 
for a non-constant potential in the t — > + limit. As an example, consider V(4>) = m y for 
which H and <f> have the following asymptotic forms at early times, 



H = — 
3t 



= 4>o 



1 - 



m 2 t 2 



+ 0(m t ) 



(2.10) 
(2.11) 
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In this case, Hip 2 ps t — > 0, so that F is essentially constant at early times. 
In this Kasner limit, the metric reduces to the following form, 

ds 2 K asner = ~<& + (Wtfidx 1 ) 2 + (dx a f (2.12) 

with p = }, $ = 0. 

This is the gravitational background in which the cosmological perturbations evolve at 
early times. The solutions for the background equations of motion suggest that the universe 
starts life with a very strong anisotropy (h — > j at early times) which is smoothed out very 
fast by the inflaton potential. The universe then enters a phase of usual isotropic inflation. 

3. Cosmological Perturbations: Small Wavelength Limit 

The computation of the spectrum for cosmological perturbations for a generic anisotropic 
background has two significant differences with the corresponding computation in the isotropic 
case: 

(1) The existence of a WKB solution for the modes of a given wavelength at asymptotically 
early times (t — > + ) is not guaranteed , since for certain backgrounds any mode may become 
super-Hubble [9] in this limit. 

(2) The £0(3) scalar and tensor perturbations are coupled for a generic wavelength at times 
t < Uso [8]. 

As commented in the previous subsection, the particular choice of the "positive branch" 
background solves (1). The positive branch metric in the Kasner limit is, in fact, a patch 
of Minkowski space-time [10] and this can be seen as follows: Let u = t sinh(Fa; 1 ) and 
v = t cosh(V r x 1 ), so that, 

ds 2 = -dt 2 + (Wtfidx 1 ) 2 + (dx a ) 2 = -dv 2 + du 2 + (dx a ) 2 (3.1) 

Therefore, this particular background admits a WKB solution for perturbations at early 
times. Note that coordinate invariants such as the curvature and the Weyl tensor are time in- 
dependent and non vanishing and hence the space is never Minkowski. Thus, when computing 
quantities that depend on derivatives of the metric one should be careful if using (3.1). 

In [8], the perturbations in the anisotropic phase were parametrized in terms of the 
variables v, H x , H + (which reduce to the usual gauge-invariant variables in the isotropic 
limit, v becomes the Mukhanov variable and H X ,H + become the two polarizations of the 
tensor modes) with the following equations of motion: 

H x +u 2 x H x =0 (3.2) 
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where the derivative is with respect to the conformal time rj and the frequencies of the coupled 
system are given as, 



"ii = e 2 '(p?+p!-2p/3 + ...) (3.4) 
"22 = e 2 P(p\+ P l- 2p$ + •••) (3.5) 

u 2 2 = ( ) (-ZW + \ (3 . 6) 

The pi are the physical momenta, p\ = kie~ 2p ,p a = k a e~ 2/3 . The ellipsis in the above equa- 
tions indicates terms subleading in the limit t — > + . To the leading order at asymptotically 
early times, e~ 2p ,uj 2 2 e~ 2p ~ 1/i 2 while uj\ 2 e~ 2p ps i 2 and as a result, the mixing terms 
can be neglected. Therefore, in this limit, the scalar and tensor perturbations decouple as in 
the isotropic case. This feature is not surprising since, as seen earlier, the metric is that of a 
flat space-time. 

In this work, we focus on the fluctuations of the scalar mode (curvature perturbation)- 
the computation for tensor perturbation can be done similarly. In the isotropic case, the 
equation of motion of the Mukhanov variable v is identical to that of a scalar field evolving 
in the same background, as long as the slow-roll conditions are obeyed (as a result of which 
^- « Therefore, in all quantities of interest, one can substitute curvature perturbation 
by a solution of the scalar field equation of motion, up to a well-defined normalization. In the 
anisotropic case, the scalar mode is given as v = exp2/3[<5</> + . 2 ^ 2 — ^ip] in terms of the 
scalar field and metric fluctuation [8]. Since v approaches the usual Mukhanov variable in 
the isotropic limit, the above argument is true for the scalar mode evolving in the anisotropic 
background in the limit t 3> ti so when the universe enters a late-time de Sitter phase. 
At asymptotically early times, v reduces to purely a fluctuation in the inflaton field in the 
Kasner background. In the limit t — > + , — ~ — — or ~ t 2 , so that v ~ 56. 
Thus, both at early and late times, the curvature perturbation can be understood as a scalar 
field evolving in the background given by the metric (2.2). Therefore, if one can find a WKB 
solution at early times, one can obtain an approximate classical solution at late times by 
a standard matching procedure at some intermediate time. This approximate classical so- 
lution will specify the particular vacuum state of the curvature perturbation field and can 
then be used to compute the late-time correlation functions. It is important to note that, in 
this scheme, the entire information of anisotropy is encoded in the vacuum state of the theory. 

Therefore, we consider a scalar field propagating in the background metric (2.2). It turns 
out that (as discussed in Appendix A) there exists a WKB solution for the scalar field in 

the Kasner regime for k 3> Hj, where Hj = and V = V(cf>)]t^o. As explained in the 
appendix, the WKB solution amounts to imposing a particular initial condition on the scalar 
modes at early times. In this work, we will only be concerned with the non-planar regime, 
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viz. k\ /c2 ~ ks, where the WKB condition is always satisfied at early times. Note that the 
condition k » Hj is equivalent to the condition of the modes being deep inside the horizon at 
early times, ^ S> 1 with e p ps y/V t and p = \. For the observed wavelength scales, between 
(1 - 10 4 ) Mpc, it can be shown that [10] 

jv-64 f T R \ A0 16 GeV \ 2 fe^ ._ 55 / T R \ / 10 16 GeV \ 2 

VlO^GevA ^ 1/4 / Hi \10 u GeVj\ J { ' 

where N is the number of e-foldings and Tr is the reheating temperature. For the range of 
observable scales satisfying the WKB condition is not hard. The visibility of the anisotropy 
at these scales, however, is not guaranteed. 



Table 1: Relations between different time scales 



Time 


Definition 


Relations 


e <*(t) 


tiso 


p(tiso) ~ PiUso) ~ tt(t iso ) 




0(10°) 




p a{U) -Ik 




> O(10 2 ) 


te 


„ait e ) — k 
e - H, 


t e ~ 2t* 


> O(10 4 ) 



At asymptotically late times [10], there exists a general solution whose precise form can 
be determined by matching the early time WKB solution at some intermediate time i*. 
The solution for the scalar field at late times can then be written as, 



= A + (k) ( f> + (ri) + A-(k) ( f>_(r l ) 



(3.8) 



where (f)±(rj) = (1 ^ikrj) exp (±ikrj), with r\ being the usual conformal time as defined in a de 
Sitter universe. 

The coefficients A + and A- are given as, 



2^2Hj . 



(2 



+ 2ie 



1 + 0(e 4 ) 



exp 



A- 



ie 3 



^jr and r 



2^/2^7 LV3 

y\kj+m 

k 



e 3 + 0(e 4 



exp 



(3.9) 



where e = 

The WKB approximation is valid when £<1 and we have retained terms up to order e 3 , 
which is the minimal order at which any signature of anisotropy appears. 
From the above solution, it follows that, 



2 

?7-S-0 



\A + + A.\ 



Hi 

2k 3 



1 + 



— r 



3/2 



cos 2\ 



k 

H~j 



(3.10) 
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Background anisotropy, therefore, implies that the ground state for the curvature perturbation 
field at late times is not given by the usual Bunch-Davies vacuum but an excited state built 
on the Bunch-Davies vacuum. 

Prom this, the late time correlation and power spectrum for the curvature perturbation 
can be easily derived. We need to replace the parameter Hi by the Hubble parameter of the 
universe at horizon exit p e and introduce the overall normalization factor of 4 also evaluated 
at horizon exit. 

Therefore, the two-point correlation function for the curvature perturbation is given as, 



Pe 2 pi 



(Ck(i)Ck'(t))«(2vr) 3 <5 3 (k + k')^3^ 



1+1 77 T 



(2,)V(k + k<)S^ 



where 




p(t e ) e p ^ w k 



cos 9 = ^- = \J\ — r 2 

k 



(3.11) 
(3.12) 



The spectral index for the curvature perturbation is then given as, 

n s -l = k^- log[F(fe, cos 9)] w log[F(£;, cos 9)] 

» 2(7] - 3e) + Q - cos 2 9^j e~ pe sin (2e pe/2 ) + 0(e~ 3pe/2 ) 

(3.13) 

Note that, e~ Pe ~ e~* e /* iso . One can easily estimate the magnitude of the correction term 
arising purely due to early-time anisotropy. In the appendix, we show that the time t* 
(time at which we match the WKB solution with the de Sitter solution) obeys t* > ti so , 
such that e p (**) 3> 1. In addition, we show that t e ps 2t*. Therefore, for t* = 5^ so , 
for example, we have t e 



= lOUso, which implies that the correction term is of the order 
10~ 5 , while the slow-roll terms are of the order 10 -2 . In this case, 
we are looking at a regime of momenta where ^- « 10 5 or, k ps 10 _1 M p (H e ~ 10" 6 M P ). 
This computation leads us to conclude that because our analytical results are valid for modes 
that exit the horizon well after the universe has isotropized, the effect of the anisotropy is 
severely suppressed in the two point function. It is important to note, however, that there 
is a corner of parameter space (3.7) (for which this analytical calculation is valid) such that 
the contribution of anisotropy to n s is large enough. For = 2ij so , we have t e = 4fj so which 
gives a correction term of the order of e~* e// * iso = e -4 ~ 1/50, comparable to the slow-roll 
terms. In this case, ps 50, such that the WKB condition is still obeyed. We will, however, 
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be interested in the momentum regime where the contribution of anisotropy to the two-point 
function is negligible and investigate its possible observable signature in the three-point func- 
tion. 

It is interesting to compare the two-point function obtained above with the ACW parametriza- 
tion [21] of the power-spectrum in a generic model of inflation with broken rotational invari- 
ance. In [21], the power spectrum was parametrized as 

P(k) = P(k) (l+g(k)(k.n) 2 ) (3.14) 

where P(k)o is the usual nearly scale-invariant contribution while n is an unit vector in a 
direction which breaks the rotational invariance. The power spectrum we have obtained is a 
slightly general form of (3.14) and admits the parametrization, 

P(k) = P(fc) (l + f(k) + g(k)(k.n) 2 ) (3.15) 

where f(k) = — | (jj^J cos (^\Z^) anC ^ = (^) COS (^\/?)' ^^ rec ^ on ^' m 
this case, can be identified with the x-axis (the scale factor along which differs from that in 
the axially symmetric orthogonal space), so that cos 9 = k.n. 

In the next section, we compute the three point correlation function for curvature pertur- 
bations in the ground state described above and analyze the effect of background anisotropy 
on non-Gaussianity of the spectrum at this level. 



4. Computation of the 3-point function 

For computing the 3-point function for the curvature perturbation, following [2], we consider 
a local (in time) non- linear field redefinition of (: 



jLr 2 + f-c 2 + 



C = C + -hC + £s£ + T^ d ~ 2 (Ccd%) (4.1) 



Evidently, this redefinition does not change the quadratic action which implies that ( c and 
( have the same equation of motion and hence the same classical solution given by equation 
(3.1). Also, since local redefinitions do not yield any enhancement of the 3-point function, it 
is sufficient to compute the correlation function in terms of the redefined field ( c . In terms of 
Cc, the leading term (in slow-roll parameter) in the interaction Hamiltonian will be given as, 

n I = -jd 3 xdr ] e^(^j P Cd~X (4.2) 

where the prime denotes derivative w.r.t. the conformal time r/ (defined in the de Sitter 
phase), and the partial indicated space derivatives. 
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We can now use the "in-in" formalism to compute the tree-level contributions to the 3-point 
function. Since there is only one kind of interaction vertex, there are only two distinct Feyn- 
man diagrams at the tree-level, viz. one with a "right" vertex and the other one with a 
"left"( recall (Q(t)) m _ m = ([Texp(i j£ Hj^dt^Q 1 (t)[T exp(-i j£ H T (t)dt)}), where T and 
T denotes the time-ordered and the anti-time-ordered product of operators. One needs to 
distinguish between vertices arising out of the time-ordered product from those coming from 
the anti-time-ordered product and we refer to them as "right" and "left" vertices respectively). 

Therefore, using the usual Feynman rules in the momentum space, the three-point cor- 
relation function at a conformal time rj is given as, 

(Cc(ki,r ? )Cc(k 2 ,r ? )Cc(k 3 ,r ? )) i?/L « ^(kj + k 2 + k 3 )A i? / L (k 1 , k 2 , k 3 , rj) (4.3) 

where A R and A L are given as, 

A R (k 1 ,k 2 M,v)=i^drj'e 3 ^ (i^j p (j2j^J jlfy'G^M) (4.4) 

A L (ki,k 2 ,ks,f7) = (A R (k 1 ,k 2 ,k 3 ,v)T (4.5) 

where G^rj,^') = Cd(ki> ??)Q(ki> rf), with ( c i(ki,rj) being the classical solution for curvature 
perturbation. A^ki, k 2 , k 3 , 77) denotes the contribution to the 3-point function correspond- 
ing to the tree-level Feynman diagram with a "right" vertex while A R (ki, k 2 , k 3 , rj) denotes 
the contribution corresponding to the tree-level Feynman diagram with a "left" vertex. The 
final result for the 3-point correlation function of the curvature perturbation is given as, 

<Cc(ki^)Cc(k2^)Cc(ks,^)) ~ <5 (3) (ki + k 2 + k 3 )[A i? (k 1 , k 2 , k 3 , rj) + A L (k 1; k 2 , k 3 , rj)] (4.6) 

In the definition of A R (k.\, k 2 , k 3 , 77), one needs to make a choice of rjo - which in standard 
inflationary scenario is taken to be —00. However, we will take rjo to be of the order of the 
isotropization time-scale when the universe has entered an essentially de Sitter phase, thus 
rjo is near the onset of inflation and parallels the choice made in [5] . In the computation 
of late-time correlation functions, the parameters p and (j) can therefore be assigned their 
respective de-Sitter values, which remain nearly constant during inflation. 
Since we are interested in the late-time correlation functions, we can set rj = 0, so that, 

G ki (r7 = 0,r/) = ^(|A + | 2 (l + iW)e- ifc ^ + |A_| 2 (l-iW)e ifci ^+ 

A+A*_(l - ikiJj'yW + A_A\{1 + ^r/)e~* M ') (4.7) 

bfGjfr = (W) = -^t^ [(\A + \ 2 + A_A\)e-^' + (|A_| 2 + A + A*_)e ik "'] (4.8) 
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In the last equation we have used n' ~ — - — 1 , n , which is valid in the de Sitter phase of 
expansion. Since we have chosen tjq ?7j so , this is a good approximation for rjo < r/ < 0. 

Now, plugging the above expression for d^G^Xr] = 0,7]') in equation (4.4), we have, 



^(k 1 ,k 2 ,k 3 )=-,5K^L f dr} > £ n^w^) 

Ac) e ±1 j =1 



(4.9) 



where the sum extends over all 8 possible linear combinations and F^ i= _i(/cj) = |^4+| 2 + 

^_^ + * and F c . = i(fci) = |^_| 2 + A\Ai*. 

Therefore, on completing the rj' integration, we have, 

A R {^MM) = -^§^- £ (U F ^))^(l-e^^) (4.10) 

Now from the expression for A±, one finds, 

F^ih) = jj (-1 + cos 2 ^ el exp (4.11) 

! + ^ f-^ + cos 2 ^ )efexp( - ) (4.12) 



The leading order term in £j is identical to what one gets in a standard computation of the 
bi-spectrum using the BD vacuum. The subleading term (« e?) carries the signature of back- 
ground anisotropy. The leading term in A R arising from the primordial anisotropy (of the 
order e 9 ) is given by the configuration £i = 1,£ 2 = — 1,^3 = —1 and its permutations. 
As is evident from equation (4.10), the bispectrum can be enhanced (in contrast to the stan- 
dard case [2]) if the denominator Yli&ki vanishes. However, the expression does not blow up 
since the exponential factor in equation (4.10) also goes to zero ( \. k . (1— exp Cih)vo) ~ 
-irjo, in the limit J2i^i -> 0). 

For the aforementioned choice of the £j,i.e. k\ = k 2 + k 3 we have the following contribution 
to the bi-spectrum, 



A(Cc(ki,0)C c (k 2 ,0)Cc(k3,0)> = (C c (ki,0)Cc(k2,0)C c (k3,0)) - (Cc(ki, 0)C c (k 2 , 0)C c (k 3 , 0)) /sotropic 

i=l 



= i <3)( Ekj) ?^-l 5 n £ «[(-i + co 8 ^ 1 ) £ ? S m 



£ 1. 

(4.13) 



with | cos B\ | = | cos 2 | = | cos 6*3 1 = cos 0. 

The above expression can now be used to estimate the contribution of background anisotropy 
to the parameter Jnl- First, let us consider the "flattened triangle" limit, where /c 2 « k 3 « 
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fci/2 ~ k. 
Naively, 

A/A _^»4( f )' (t , 0) (4 + co ^) sm H (4 . 14) 

where (Cc(ki, r?)Cc(k 2 , r?)Cc(k 3 , r?)) = (ki + k 2 + k 3 )^(k 1; k 2 , k 3 , r?) and (Ck(*)Ck'(*)> = 
(27r) 3 (5 3 (k + k')P(k). However, as pointed out in [5], a factor of \krjo\ is lost when one 
computes the l-space bi-spectrum. Therefore, the final contribution to Jnl becomes, 



Af NL = e^j ^__ + cos^J (4.15) 

which, in addition to the standard slow-roll factor, is suppressed by powers of |, leading to 
an extremely small change in the value for /nl- 

Now, consider a "squeezed" limit of this configuration: fc 3 <C & 2 ps k\ ~ A;. In this case, 

A^(k,k,k 3 ) P{p\K. ,f 1 2 ^ ,, 1fi , 

A/iVi ~ p(k)p ( k 3 ) " ^ UJ ( %) v 3 + cos V (416) 

leading to the same final contribution to Jnl as derived in the "flattened triangle" limit. 
Equation (4.10) can be studied in another interesting limit , namely k 3 <C k\ ps fc 2 ~ k but 
l^i — ^2 1 / ^3 - usually known as the "squeezed triangle" limit. Note that the denominator 
Si Ciki = —ks, where | £13770 1 S> 1 ensuring the mode is sub- horizon around the time when the 
universe isotropizes. 

The resultant contribution to the bi-spectrum is given as, 

AA( kl ,k 2 ,k 3 ) « ^-f^ ^n^ [(-^ + COs2 ^) ^ + (-i + cos^ 2 ) 4^ 

' +c.c. (4.17) 



£3 

Setting /ci « fc 2 ~ in the above equation, we have, 
AA(k 1 .k 2 ,k 3 ) « £ C| )' (^)" 2 (n ^ftW»-2/S) ' H -? + ^ + *"»>' (418) 



/,2p3 



Therefore, 



3 



AA(k,k,k 3 ) _ (py k 

p(k)p(k3) ~ C ^j ^ (4.19) 

which shows that the enhancement for the "squeezed triangle" limit is greater compared to 
the "flattened triangle" by a factor of k/k 3 . The result is in agreement with the findings 
of [16], where the bi-spectrum in the "squeezed" limit was shown to be greater than that 
obtained in the "flattened" limit by this precise factor of k/k$, for a generic (i.e. non-Bunch 
Davies) choice of the vacuum state. It was also argued in [16] that (k/ks) 

max ~ 200, given 
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the range of scales measurable in the CMB. For e ~ 10 , p ~ Hi « 10 M p and comoving 
wave-number k ps 10 _4 M p , we have |A/jvl| ~ 10~ 3 , which shows that the "squeezed triangle" 
enhancement is too small to produce any observable non-gaussianity. 

As explained in [5], in order to have a larger Jnl, one needs to include terms involving 
higher derivatives in the scalar field, which are compatible with the slow-roll conditions. As 
explained in [7], the lowest dimensional operator that can be included is a dimension 8 oper- 
ator , where M is the cut-off scale for the effective field theory for inflation. As a next 
step, we compute the contribution of the background anisotropy once the higher derivative 
interaction terms are also included. The computation roughly follows the treatment in [5]. 
The single-field theory,analyzed above, should now be modified by adding a term Cj = 
- v /^7g^ ¥ (V</>) 4 . On expanding the lagrangian around the classical solution to third order 
in perturbation, we have the following interaction Hamiltonian for the curvature perturba- 
tion: 

Ui = - J d 3 xa(r,)^j^C'(C' 2 ~ (^C) 2 ) (4.20) 

As explained in [5], the three-point correlation function at the tree level (from the vertices 
above), computed for an excited state, is enhanced for a flattened triangle configuration, as 
we had observed in the case of a single-field theory without higher derivative interaction. In 
addition, we shall have an enhancement for a "squeezed triangle" configuration as well. 

From (4.20), one can directly compute the 3-point correlation function, 



A^kx.ka.ks) = i [° dve p ^^[Hd v G kl (0, V ) x (3!) 

Jrio P i=i 



+ ((k 1 £ 2 )G kl (0,r ] )G k2 (0,r ] )d ri G k3 (0, V )+ perms) x (2!)] + c.c. (4.21) 

where the factors of 3! and 2! are the respective combinatorial factors for the two vertices. 
As before, the leading terms due to anisotropy will appear in the integrand as coefficients of 
e (Ei5» fc l )»? ) w ith £i = 1,£ 2 = —1,^3 = — 1 an d its permutations. It is, therefore, sufficient to 
isolate the contribution proportional to exp ^ f (with £jS specified above) in the integrand 
for computing enhancements in the bi-spectrum in the flattened as well as the squeezed limit. 
From the definition of Yli=i 9f]G ki (Q,rj), we have, 

Y[d v G k .(0,v) = (£)V jff [ e ^i- fc 2- fc 3)r; g 3 ( _ 1/3 + cog2 g^-a/ei + perms ^ ( 4 . 22 ) 



k 3 r]- 



16p 3 



e i( kl -k 2 -k 3 ) v ^ _ ^ + ik2r] + klk27] ^ £ 3 



(-1/3 + cos 2 6>i)e 2i/£l +perms. 



(4.23) 
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Plugging these back into (4.21), we obtain 

A R = i f dn -^- A in ^f ] (|) 6 [e^- k ^el(-l/3 + cos 2 e^e'^A^MM, n) 

+perms.] + c.c. (4.24) 
where we have used e p r] = 1/p. The explicit form for the function »4(ki, k2, k3, rj) is given as, 
^(ki,k 2 ,k3,»7) = ^(JJfc?) + {kj -k\- k 2 )(l - ik lV + ik 2 v + hk 2V 2 )k 2 

i 

+ (k 2 -k\- k\){\ + ik 2 r] + ik 3 r] - k 2 k 3 i] 2 )kf 

+ (k 2 -k\ - fc|)(l - ikir) + ik 3 rj + kik 3 r) 2 )k\ (4.25) 

From the above general expression, we can now proceed to compute Jnl in the limits of 
interest. 

Consider first the flattened triangle limit: k 2 = k 3 = k = fci/2. The function -4(ki, k 2 , k3, r?) 
reduces to, 

^4(ki, k 2 , k 3 , 77) = 24&V + *fc^(4fe 4 + 16fc 4 + 4fc 4 ) + (-4A; 4 + 8fc 4 - 4A; 4 ) + k 2 r, 2 {-8k 2 - 8k 2 - 8k 2 ) 
= 24ik 5 r] (4.26) 

Integrating out rj, we have, 

Ar « ^ I ^ £ 3 r ?0 2 cos2/£(-l/3 + cos^) (4.27) 

Therefore, the ratio of the enhancement of the bi-spectrum in the flattened triangle 
configuration to the bi-spectrum computed in the standard Bunch-Davis vacuum is given as, 



^(C(ki,r]o)C(k 2 ,Vo)C(k3,Vo) / / \ 

e 3 _ +C os 2 )\kno\ 2 (4.28) 



C(4%)C(4%)C(4%)) Dn V 3 

where A (C(ki, 77o)C(l £ 2 , ?7o)Cc(k3, r?o)) denotes the leading contribution of background anisotropy 
to the bi-spectrum as before. 

However, as in the previous case, one factor of \krjo\ drops out from the ^-space bi-spectrum 
and as a result, the /nl parameter for these configurations changes by, 

IA/atlI « f^Nol (4.29) 

where e = |&f?o| = a (ri )H(r)o) — TTi an< ^ * s re ^ a ^ e< ^ to the slow-roll parameter e as 

0e ~ Pe 2 M 2 e (here we have replaced the parameters <j) and p by their values at horizon exit, 
assuming that these remain essentially constant in the de Sitter phase). 
This implies that, 

. . „ . Hj Mp e e 3 

IA/atlI - M P 3 (4.30) 
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To estimate the range of values for /ml, we take Hi « H e « « H(j]q) 10 6 M P , 
where 77n > 77j so . We define a number n such that n = t^ 2 -. 

Now, for the range of momenta we are interested in, we require fce~ a<T *) < M, where t* cor- 
responds to the time at which the early-time WKB solution is matched with the late-time de 
Sitter solution (see appendix) and is defined as e a(T *) = a/ 77-- Combining the two relations, 



we clearly have k 1 / 2 < -jjfa. Taking H* « Hi, we find e = sj S± > 

Therefore, plugging in the above in equation (4.30), we have, 

^ Hi fM p \ 3 (Hi\ 3 
\Af NL \ « — — e 



M p \M) e \M) 



M„7 V M 



Now, substituting the values for the slow-roll parameter e ~ 10 and the scale of inflation 
Hi « 10~ 6 M p , the above equation implies that a large /nl will require M ~ 10~ 4 — 10 _5 M p . 
For example, Jnl ~ 100 if we choose M = 2.5 x 10~ 5 M P . 

Next, we consider a "squeezed triangle" limit which was not addressed in [5]: k 3 <C fci ps 
^2 ~ A;, with fc] ^fe + ^3. In this limit, the function .A(ki, k2, k3, 77) reduces to, 

^(ki, k 2 , k 3 , 77) = 6k A k 2 3V 2 + (-2k 2 )(l + fcV )*§ + (-fcf)(l + ikr) - kk 3V 2 )k 2 
+ (-k1)(l - ikr] + kk 3 i] 2 )k 2 

= -4k 2 k 2 + 4fc 4 r ? 2 fcl (4.32) 

(4.33) 

Now, 



r0 rO 

/ d^^-^-^i-Aekl+A^klr) 2 ) « / dr/e- ife3, '(-4fc 2 fc|+4fc 4 fc|r ? 2 
J 170 ■/ 170 



„3 / fc M .,,3 / M - A1 Ji ( k 



(4.34) 

since k ^> k 3 . In evaluating the integral we have assumed that I&377I 3> 1 for r\ £ [770 , 0] . 
From (4.24), one obtains the following bispectrum, 

A R « ^i^4 e3 cos 2/e(-2/3 + cos 2 ^ + cos 2 2 ) (4.35) 
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where e = \/|. 



The corresponding /jvl can be read off as follows: 



•Je 3 / k \ f Hj \ ° f M p \ 1 k 



Evidently, 



I a rsqueezed\ 



. =|%?o|»l (4.37) 

'7VL 

which shows that enhancement in the bi-spectrum for flattened triangle limit can be more 
important than squeezed triangle limit, once we include higher-derivative terms in the inter- 
action Hamiltonian. 

To estimate the range of f^ eezed , we need to ascertain the maximum value of ^. Recall 
that \kr)o\ < Therefore, taking M ~ 10 _4 M P , we may set |fc?7o| max ~ 100. Now, since 
|&3?7o| S> 1, one can also set \k3ri0\min ~ 10. This naturally fixes the ratio (j^) ma x ~ 10 - 
one order lower than the estimated value in [16]. In the squeezed limit, therefore, the scale 
M cannot be lowered much below 10~ 4 M P , as it pushes the ratio (£-) m ax closer to 0(1). 
Plugging these values in (4.36) we have, \Af^ eezed \ « 10~ 3 , which is too small compared to 
the observable limits of non-gaussianity. 

This proves that the flattened triangle is the dominant source of non-Gaussianity for the 
curvature perturbation spectrum given an interaction Hamiltonian of the form (4.20). 

Finally, the particular value of the cut-off scale M for the effective field theory merits 
some explanation. In [14], it was argued that the scale M cannot be too small compared to 
\/2eM p if one wishes to impose a limit on the size of higher derivative terms in the effective 
action. For M = 5 x 10~ 4 M P , however, one cannot rule out the presence of higher dimensional 
operators in the effective action and their contributions to A/jVL- These contributions will 
be suppressed by factors of ^ [7] where ~ 10 -2 — 10 _1 in our case (as opposed to a 
suppression factor of 10" 5 in [14]) 



5. Conclusion 

In this note, we have presented a way to probe the signature of pre-inflationary background 
anisotropy in the spectrum of cosmological perturbations at late times for a range of high- 
momentum non-planar modes. These modes, which generically exit the horizon after the uni- 
verse has isotropized, have a nice WKB description at early times. The late time dynamics of 
these modes is characterized by an excited state built on the standard Bunch-Davies vacuum 
and this state carries the signature of the pre-inflationary anisotropy. We have computed 
correlation functions of the curvature perturbation (two-point and three-point functions) in 
this new ground state and investigated the issue of having possible observable signatures of 
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anisotropy. 



Our computation suggests that contribution of anisotropy to the spectral index could be 
appreciable even for the range of large momenta modes (for which the WKB approximation 
is valid) which exit the horizon after isotropization. However, if the number n = -M- is of the 
order of 10, any correction to the two-point function is severely suppressed in a single-field 
inflationary model with a canonical action. The effect of anisotropy will obviously be much 
larger for modes with small wavenumbers which exit the horizon before isotropization. The 
two-point function in this regime was numerically analyzed in [8]. 

The three-point function for the curvature perturbation, in this case, is enhanced in a 
flattened triangle as well as a squeezed triangle limit, although the latter leads to a larger 
Jnl ■ However, the bispectrum in either case is suppressed by factors of in addition to the 
usual suppression by the slow-roll parameter, making the contribution negligible compared 
to observable limits. 

The problem can be circumvented, to a certain extent, by including higher derivative 
operators in the action and in this work, we have studied the effect of having a dimension 
8 operator g^r(V0) 4 in the action. Computation of the resultant three-point correlation 
function shows that the bi-spectrum for the curvature perturbation is again enhanced for a 

flattened triangle configuration and a squeezed triangle configuration, with the former being 

\Aff lat \ 

the dominant source of non-Gaussianity in this case - ~ 10 for the cut-off scale 

M ~ 10~ 4 Mp. 

We have shown that the /nl for the squeezed triangle case will always be extremely 
small, since, for reasons explained at the end of the previous section, one cannot lower the 
cut-off scale much beyond 10 _4 M p . For a large f^L ~ 100, in the flattened triangle case, one 
needs to have the cut-off scale for the effective field theory to be set around M ps 10~ 5 M p , 
which is low enough for higher-dimensional operators (suppressed by factors of ^ ) to appear 
in the effective action. We observe that this result is very similar to that obtained in [5] where 
the authors studied the enhancement of the bi-spectrum for an excited state in a single-field 
theory with a dimension 8 operator. 

In our case, however, the effective ground state is characterized by an angular dependence 

and a particular momentum dependence (~ (y^f) 3 )- These are the distinctive signatures 
of background anisotropy in the enhanced non-Gaussianity for the flattened triangle and the 
squeezed triangle configurations, distinguishing it from a generic case of enhanced bi-spectrum 
for a non-BD vacuum state. 

Finally, we want to remind the reader that this analysis holds only for the non-planar modes, 
i.e. modes with k\ ~ k2,k% for which the WKB condition always holds at early times. The 
planar case and non-Gaussianities from the interplay of these two regimes are the subjects of 



- 16 - 



a work in progress. 
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7. Appendix: Scalar Field in Anisotropic Background 

In this section, we study the evolution of a single, massless scalar field minimally coupled to 
gravity in an axially symmetric anisotropic space-time, with a positive cosmological constant. 
In particular, we construct a WKB solution for the scalar field at early times and show how 
one can match it with the general solution at late times to obtain equations (3.9). The 
treatment essentially follows [10] 
We consider the following action for the scalar field, 

S = -J d^x^-g Q<T^<R0 + V) , (M 2 V = 1) (7.1) 

where the background metric is chosen to be an axially symmetric version of the Bianchi I 
metric: 

ds 2 = -dt 2 + exp(2p)(dx 1 ) 2 + exp(2/3)(dx a ) 2 , (a = 2, 3) (7.2) 
where p, j3 are known functions of time: 

p = ^lntanh 2 (^J^j sinh(3F/£) 



sinh (SHrt) 

3 -¥) 



tanh 



(7.3) 



with Hj = y -j . 

Define p = a — 2/3+, /3 = a + /3+ and a new "time" coordinate r, analogous to the conformal 
time in the isotropic limit, as, 

dr = % (7.4) 

From equation (7.3), one can derive, e 3a = e p+2 @ = smh(3Hjt) = sin fr(_3_f/ fr ) • 
It can be easily seen that as t varies from + to oo, r varies from — oo to 0~. In this time 
coordinate, the equation of motion for a mode 4>k is given by, 

^+^(t) 2 )<^ = (7.5) 
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The frequency squared is given as, 

u (r) 2 = ^-{l-^x) (7.6) 

where r 2 = k2 ^ 3 and x(t) = 1 — e 6H ' T = exp (— 6a)(-y/exp (6a) + 1 — 1). Evidently, x(t) 
varies from 1 to as r changes from — oo to 0. 
Equation (7.6) has a WKB solution: 

1 f T i 

4>wkb = —j= exp \—i \ dr Cj] (7.7) 
V 2tD J r 

where uj has to be determined from the equation, 

The WKB approximation holds as long as the WKB parameter 

dw 2 
dr 



W 3 



Hi 1 - x(i 



+ 1 < 1 (7.9) 



k (x/2) 1 /3(l_ r 2 x ( r ))l/2 ^l_ r 2 x ( r ) 

The choice of the WKB solution above is obviously equivalent to imposing a particular initial 
condition on the modes of the scalar field at early times. This can be seen directly by 
analyzing the early time behavior of the classical solution for the scalar field. Firstly note 
that oj ~ uj — > k\ in the limit r — > — oo (or t — > 0+) and in this limit r and t are related as 
r = jjtj In 3 ^ 7 * . Therefore, in the early time limit, the time-dependence of the WKB solution 
is given as follows, 

1 fci , 

4>wkb ~ -7= ex P[-^^rr lo g — ^— J ~ * 3H/ (7-10) 
Vzfci oil/ z 

Now consider the equation of motion of the scalar field at early times, 

1 k 2 
4>k + -(t>k + {kl + kl + r^)4> k = ^ (7.11) 

Define z = In so that the equation reduces to 

^W + (e 22 + ^)0^) = O (7.12) 

In the limit t — > 0+, z — > —00, so that the exponential term drops out of the above equation 
and we obtain a solution of the form, 

<f> k (z) = ^(fc)e - ^ 2 + B{k)e~Jsv z (7.13) 
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Choosing B(k) = 0, we find that (j>k(z) ~ t 3H i , confirming that the WKB solution has 
the same time-dependence at early times as expected from the classical solution subject to a 
certain initial condition. 

In fact, equation (7.11) has the general solution, 

cj> k (t) = C^H^^yJq + h^t) + C 2 (k)HW /y/w {yJk* + k*t) (7.14) 

Using H^(z) = —^(J^e- - J_ iv ( z )) and H^(z) = —^(-J^-™ + J_ iv (z)), 
we can rewrite the general solution as, 

Mt) = A(k)J iki/VW (y/k* + k*t) + B{k)J_ iki/VW {^kl + klt) (7.15) 

Therefore, the WKB solution corresponds to imposing the initial condition A(k) = and 
choosing B(k) appropriately, as ^_jfc 1 / v ^y(\/^2 + ^P) nas ^ ne same time-dependence at early 
times as our WKB solution. 



We will be interested in the large momentum regime of non-planar wavenumbers, i.e. 
ki 3> H, implying that the factor (1 — r 2 x(r)) in the denominator doesn't vanish anywhere 
(since both r and x are fractions). In the regime x ~ 1 (i.e. early times), the WKB condition 
obviously holds for any momentum. In [10], it was shown that the condition holds for high- 
momentum modes as long as 3> expa(t). Therefore, the time at which the WKB solution 
should be matched with the late-time de Sitter solution has a natural choice, r*, such that, 

e»( T *) = ./A (7.16) 
V Hj 

The above equation implicitly states that r* corresponds to late times when e a ( T *) 3> 1. In 
terms of real time, this condition implies that e a ^*^ ~ e Hlt * 3> 1, or t* > ij so ~ Also, if 
t e denotes the time of horizon exit for a given mode of wavenumber k, we have e Hjte = jj. 
This suggests a simple relation between i* and t e , viz. 

t e « 2t* (7.17) 

Returning to the problem of matching the modes - since k 3> Hj, x(r*) « exp (— 6a(r*)) ps 
0. Therefore the WKB solution around r = r* is given by expanding equation (7.7) around 
x = and then plugging in the values of x and the frequencies at t = t*. The solution can 
be expanded in powers of e = \J~^f and one needs to retain terms to the order at which the 
direction dependence first appears. It turns out that it is sufficient to retain terms up to the 
order e 3 and to this order x(t#) and the frequencies are given as, 

x{U) = 2(H I /kf 2 (l - {Hj/kf 2 ) (7.18) 
k 2 



= — 
H 1 



1 + ( F r2)( f )3/2 



(7.19) 
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ui = ut i - 2 



.fir 

Now, in the de Sitter regime the solution to equation (7.5) is given as, 

(j) k = A + + (t) + A_0_(t) 
where the modes <p± are given as, 



ik 



±(t) = 1 =F ^(-Sff/r) 1 ^ exp ±^(-3^r) 1 /3 



#7 



#7 



(7.20) 



(7.21) 



(7.22) 



Matching the de Sitter solution with the WKB solution at r = t*, we obtain the following 
equations for and ^4_, 



A^' + (t*) + ,4_0'_(t*) 



A+^+(t*) + A_0_(t*) = 1 + 



#7 




1 



2k 



2k 



1 

2uj* 



(7.23) 
(7.24) 



where we have absorbed an overall phase in the definition of and A-. Solving for A± 
from the above equations, we have, 



A + = 



IE 



2^2H} 



A- = 



IE 



2^/mi 



(2- £ 2 )+2* £ (^-l)+0(e 4 ) 
2 



exp l — 



r 2 \e a + 0(e a ) 



exp 



One can easily verify that these coefficients obey the normalization condition, 



(7.25) 
(7.26) 

(7.27) 
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